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Abstract 

The neutralino relic abundance is evaluated for a wide range of the neu- 
tralino mass, 20 GeV < m x < 1 TeV, by taking into account the full set 
of final states in the neutralino-neutralino annihilation. The analysis is per- 
formed in the Minimal SuSy Standard Model; it is not restricted by stringent 
GUT assumptions but only constrained by present experimental bounds. We 
also discuss phenomenological aspects which are employed in the companion 
paper (II. Direct Detection) where the chances for a successful search for 
dark matter neutralino are investigated. 
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1. Introduction. 

The general features of cosmological structures, as they are observed and 
understood at present, lead to the conclusion that large amount of the matter 
in our Universe is in the form of Cold Dark Matter (CDM). This circumstance 
has recently prompted new detailed investigations about the neutralino (x) , 
since this SuSy particle appears to be the most favorite candidate for CDM. 
Here we analyse one of its basic properties, the relic abundance, by extend- 
ing a previous investigation of ours [1] (which was confined to neutralino 
masses below the W-boson mass) to a much wider range of neutralino masses: 
20 GeV < m x < 1 TeV for the most general neutralino composition. Fur- 
thermore we take into account the whole set of exchange diagrams and final 
states in the x ~ X annihilation process which is the fundamental ingredient 
in the evaluation of the relic abundance. In addition, radiative corrections 
to the Higgs boson masses as well as to the relevant coupling constants are 
appropriately included in our evaluations. 

In the present investigation the theoretical framework is represented 
by the Minimal Supersymmetric Standard Model (MSSM); only one stan- 
dard GUT assumption is employed with the purpose of simplifying the phe- 
nomenological discussion. 

The analysis presented here differs from the previous ones [2-8] in at least 
one of the features mentioned above. In particular the most recent analyses 
on this subject by other authors are mainly based upon GUT schemes which 
also include supergravity; this automatically implies model-dependent rela- 
tionships between the various masses which come into play. Here we prefer 
to consider a more flexible scheme where unknown masses are not a priori 
fixed, but are only constrained by present experimental bounds. 

We also wish to emphasize that we do not restrain our attention to re- 
gions of the parameter space where the neutralino would provide by itself the 
total amount of required CDM, but we rather widely explore the parameter 
space with the aim of investigating the chances to detect the neutralino as a 
dark matter candidate by direct or indirect searches. In fact regions where 
the detection event rates are higher do not necessarily coincide with the lo- 
cations in parameter space where the relic density is larger, since stronger 
coupling of the neutralino with matter may compensate for partial deple- 
tion in the neutralino local density. Indeed actual detection of dark matter 
neutralino would be an achievement of paramount interest even if the neu- 
tralino does not exhaust our need for CDM. For these reasons, whereas in the 
present paper (hereafter called I) we discuss the neutralino relic abundance, 
in the companion paper (II) which follows this one, we analyse the problem 
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of its direct detection. The way of presenting our results in paper I is mainly 
shaped according to the needs for the applications discussed in paper II. 

2. Minimal SuSy Standard Model. 

Our theoretical framework is the MSSM, with the standard definition 
of neutralino as the lowest-mass linear combination of photino, zino and 
higgsinos, 



where 7 and Z are linear combinations of the U(l) and SU(2) neutral gaug- 
inos, B and W3, 



9w being the Weinberg angle. As usual, x is assumed to be the lightest super- 
symmetric particle (LSP) and then stable if R-parity is conserved. Its mass 
as well as its composition depend on the parameters: Mi, M 2 (masses of B 
and of W3, respectively), [i (Higgs mixing parameter) and tan/? = v u /vd (v u 
and Vd being the v.e.v.'s which give masses to up-type and down-type quarks). 
A standard procedure in current literature is to embed the MSSM in GUT, so 
that a relationship between Mi and M 2 follows: Mi = 5/3 tan 2 6\y — 0.5 M 2 . 
Relaxing this assumption may modify some neutralino properties in a sig- 
nificant way; this has been discussed by some authors in the restricted case 
m x < mw [9,10]. 

In this paper we will report our results under the assumption that the 
GUT-induced relation between Mi and M 2 holds, as this is a most natural 
hypothesis. Extensions of the present analysis due to the relaxation of this 
relationship will be presented in a forthcoming paper [11]. 

As mentioned above, in the present analysis also high values of m x are 
considered (up to 1 TeV). Thus, in order to deal with a parameter space large 
enough to contain various compositions for x at any value of m x , we have 
taken wide ranges for M 2 and fj, : 20 GeV < M 2 < 6 TeV, 20 GeV < < 
3 TeV. The scatter plots which are presented in section 3 are obtained by 
varying M 2 and \i in these ranges. As far as tan (3 is concerned, in order to 
cover a wide range for it we will representatively choose the following values: 
tan/? = 2, tan/3 = 8 (or rather the range 6 < tan/3 < 10 in a number of 



X = ai7 + a 2 Z + a 3 H® + a 4 -ff 2 



(2.1) 



7 = cos 9\yB + sm$wW3, 
Z = — sin 9wB + cos 9wW 3 , 



(2.2) 
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extended scatter plots) and tan/3 = 20. It is worth noticing that the ranges 
chosen for the parameters automatically disregard regions of the parameters 
space which have been excluded by LEP. 

To illustrate the general features of the neutralino we give in Fig.l a 
representation of its mass and composition in the Mi — \i plane at tan/3 = 8. 
This figure clearly shows that the neutralino composition tends to be very 
pure (either pure higgsino or pure gaugino, depending on whether Mi > 2\n\ 
or Mi < 2|/i|) as M 2 and p increase. This property, already discussed in 
previous works, is due to the fact that in the neutralino 4x4 mass matrix 
the higgsino sector asymptotically decouples from the gaugino sector as Mi, 
\p\ >> Mz- Before this asymptotic regime turns on, and so for M 2 , 2|/z| < a 
few hundreds GeV (neutralino masses < 100 GeV), mixed higgsino-gaugino 
configurations are important. This property plays a crucial role in some of 
the conclusions we will draw about the neutralino relic abundance and the 
detection event rates. 

Effective interaction properties of the neutralino depend sensitively on 
the three physical neutral Higgs bosons of the MSSM: the two CP-even 
bosons: h,H (of masses rrih, ran with ran > rah) and the CP-odd one: 
A (of mass m^). Radiatively corrected relationships between these masses 
are used here. Thus, taking rah as an independent parameter, ra A and ran 
are functions of m^, tan/3, rat (mass of the top quark) and fa (mass of the 
top scalar partners, taken as degenerate). To represent our results we will 
take two values for m^: rah = 50 GeV (which entails ra A ~ 50 GeV) and 
ra h = 80 GeV (this implies ra A = 83 GeV for tan/3 = 8 and ra A = 170 GeV 
for tan/3 = 2). The top quark mass has been taken at the value ra t = 150 
GeV. The sfermion masses will be discussed in the next section. 

3. Relic Abundance. 

The neutralino relic abundance £l x h 2 (O is the density parameter of the 
Universe: O = p/p c , p c = 1-88 x 10~ 29 h 2 g cm -3 , 0.4 < h < 1) has been 
evaluated employing the usual formula [12] 



^, 1 3x 10 -(|) S (^)V(;^). (3-D 

Here T 7 is the present temperature of the microwave background, T x /T 7 is 
the reheating factor of the photon temperature as compared to the neutralino 



-4- 



temperature, Xf = Tf/m x ~ 1/20 where Tf is the neutralino freeze-out 
temperature. The ^/-dependent expression in the denominator of Eq.(3.1) 
represents the integration from Tf down to the present temperature of the 
thermally averaged quantity < a ann v >= a + bx where a ann is the x ~ X 
annihilation cross-section, v the relative velocity and x = T/m x . In the 
evaluation of the annihilation cross-section we have considered the whole set 
of final states: 1) fermion-antifermion pair, 2) pair of neutral Higgs bosons, 
3) pair of charged Higgs bosons, 4) one Higgs boson-one gauge boson, 5) 
pair of gauge bosons (W + W~ , ZZ). For the final state 1), the following 
diagrams have been considered: Higgs-exchange diagrams and Z-exchange 
diagram in the s-channel, /-exchange diagrams in the t-channel. For the 
final states 2-5) we have taken Higgs-exchange and Z-exchange diagrams in 
the s-channel, and either neutralinos (the full set of the four mass eigenstates) 
or chargino exchange in the t-channel, depending on the electric charges 
of the final particles. As for the sfermion masses, we have considered two 
extreme cases: one with mass values as low as possible, compatibly with the 
present experimental bound and with the assumption that x is the LSP, and 
a second one where all sfermions are very massive. The lower limit that we 
have conservatively used for the sfermion masses is the LEP bound of 45 
GeV, since the CDF limit [13] on the squark masses does not appear to be 
consistent with a massive neutralino [14] . 

Thus we have considered the two cases: 1) mj = 1.2 m x , when m x > 45 
GeV; mj = 45 GeV otherwise, except for the mass of the top scalar partner 
(the only one relevant to radiative corrections) which has been taken fa = 3 
TeV ; 2) mj = 3 TeV for all sfermions. 

Let us now turn to the presentation of our results. In Fig.s 2-3 we report, 
in the Mi — plane, the dominance of various final states at tan/? = 8, 
rah = 50 GeV and for the two representative values for mj. Dominance of 
a particular final state means that this channel weighs for at least a factor 
of 5 over the other states in the quantity < a ann v >. Fig. 2, which refers 
to small values of raj, shows that in large regions of the M2, fx plane the 
// and the two gauge bosons final states dominate. However it has to be 
stressed that for \x > in a significant region of the parameter space, where 
mixed gaugino-higgsino configurations occur, there is dominance of one Higgs 
boson-one gauge boson final state. This feature was absent in Ref. [3], 
where this contribution was estimated to be subdominant. The one Higgs 
boson-one gauge boson final states are even more important when all the 
raj masses are set at 3 TeV (see Fig. 3); in fact, a drastic increase in mj has 

the obvious consequence of a suppression of the /-exchange amplitude, and 
this significantly depresses the // final states. On the contrary, dominance 
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of the two Higgs boson final states is limited to a few isolated points in the 
parameter space, as is illustrated in Fig.s 2-3. These general features remain 
almost unaltered if we move from tan/? = 8 to smaller values, tan/? ~ 2. 

Let us now discuss the neutralino relic abundance. In Fig. 4 we report 
our results in the form of scatter plots obtained by varying M2, \i in the 
parameter space previously defined and by varying tan/3 in the range 6 < 
tan/3 < 10. O x /i 2 versus m x is shown for three different types of neutralino 
compositions: higgsino dominance (dominance here means 90% or more), 
gaugino dominance, maximal higgsino-gaugino mixing (i.e., 0.45 < a 2 + a| < 
0.55). In Fig. 4a, (higgsino dominance) some characteristic features are very 
clearly displayed. A number of pronounced dips (and of sharp falls off) in 
£l x h 2 reflect the presence of poles (and the opening of new thresholds) in the 
annihilation cross section. In sequence we have: at m x ~ 25 GeV the h and 
the A poles, at m x ~ 45 GeV the Z pole, at m x ~ 90 GeV the threshold 
for the x ~ X annihilation into channels W + W~ and ZZ. In the case of 
gaugino dominance (see Fig. 4b) sfermion-exchange amplitudes provide large 
contributions to the annihilation cross section, with the effect of depressing 
the neutralino relic abundance as compared to the case of higgsino dominance 
for m x ~ 90 GeV. At higher m x values the gaugino-dominated compositions 
give a larger Q x h 2 , since here the // final state is dominant, but somewhat 
hampered by the running values of the / mass: mj = 1.2 m x . Compositions 
with large mixings (see Fig.4c-d) entail rather low values of Q x h 2 , due to 
the substantial contribution to the annihilation provided by Higgs-exchange 
and /-exchange. Nevertheless, these neutralino configurations contribute 
significantly to the event rates for direct neutralino search (as discussed in 
paper II). The previous discussion should make the features of Fig. 5 quite 
transparent. In fact here high values for mj force the neutralino relic density 

of the gaugino-dominated configurations to be large, by inhibiting the /- 
exchange amplitude. This same mechanism is the reason for the depletion 
in the // final state dominance that we notice in the plot of Fig. 3 when we 
compare it with the one in Fig. 2. 

Fig.s 6-7 show the neutralino relic abundance when tan/3 is small: 
(tan/3 = 2, for definiteness) and nth = 80 GeV. Higgsino-dominated con- 
figurations at small m x (below the thresholds for W + W~ and ZZ) display 
large values of fi x /i 2 , since now, smaller values of tan/3 and larger values 
of Higgs boson masses, both have the effect of suppressing the annihilation 
channels with Higgs exchanges. As for the gaugino-dominated compositions 
it is worth noticing that in the case of large mj (see Fig. 7b) Q x h 2 displays 
a behaviour which is rather common in the context of some supergravity 
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inspired models. In fact in these schemes it frequently occurs that theoreti- 
cal and phenomenological constraints restrict tan (3 to very small values and 
sfermion masses to high values, with the consequence that the cosmological 
requirement Q x h 2 < 1 can only be met at the Higgs-poles or at the Z-pole. 
Consequently, particular care has to be taken in the evaluation of the relic 
abundance [15] in these models. In our kind of analysis, fine-tuning of m x 
with the masses of the Higgs bosons or of the gauge bosons is not required 
and would then appear rather accidental. 

A final scatter plot for O x /i 2 is shown in Fig. 8 for tan/? = 20 and m>h = 
50 GeV. 

A word of warning is required about the effects due to the possible 
occurrence of an approximate degeneracy (within ~ 15%) between the neu- 
tralino and some other SuSy particle. When this happens the neutralino 
decoupling mechanism is enhanced due to the annihilation process involving 
the neutralino with the other SuSy particle which is close to it in mass (this 
process is usually denoted as coannihilation in the literature) [15,16]. Ef- 
fects on O x /i 2 due to coannihilation may be large (one order of magnitude or 
more, depending on the nature of the coannihilating particle and on other de- 
tails of the theoretical scheme) . Apart from the peculiar and accidental case 
when for instance a sfermion and the x would almost have the same mass, 
a natural case of approximate degeneracy occurs in the neutralino-chargino 
sector. However, this happens in regions of the parameter space of higgsino 
dominance. In paper II it is shown that chances of detecting dark matter 
neutralinos rely essentially on mixed or gaugino compositions for neutralinos. 
Thus coannihilation does not significantly affect the evaluations of the event 
rates presented in paper II. 

In conclusion our analysis confirms that the neutralino, even at mass 
values higher than the W-mass, may satisfy the attributes required for a 
good candidate for CDM. It is also clear that the theoretical evaluations for 
neutralino relic abundance only suffer from the lack of information about 
some of the particles that naturally come into play, such as the Higgs bosons 
and any SuSy object. Only new experimental inputs can help theory in 
sharpening its predictions for neutralino dark matter. 

* * * 

This work was supported in part by Research Funds of the Ministero 
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Figure Captions 

Figure 1. Isomass curves and composition lines for neutralino in the 
M 2 — fi plane for tan/3 = 8. Dashed lines are lines of constant x mass (m x 
= 30 GeV, 100 GeV, 300 GeV and 1 TeV). Solid lines refer to constant 
gaugino fraction f g in the neutralino composition (f g = a\ + a 2 ): f g = 0.99, 
0.9, 0.5, 0.1 and 0.01. 

Figure 2. Final states dominance regions in < a ann v >, for tan/3 = 8 
and rrih = 50 GeV. Sfermion masses are given by: mj = 45 GeV, when 
m x < 45 GeV; mj = 1.2 m x otherwise (except for the SuSy partners of the 
top quark whose common mass is set at m = 3 TeV) . Dominance region for 
a particular channel is defined as the region where that channel dominates 
over the other ones by a factor of five at least. Different regions are marked 
as follows: heavy dots for // final state, horizontal lines for gauge boson 
pair final state, diamonds for mixed Higgs boson - gauge boson final state, 
squares for Higgs boson pairs final state. In regions marked with light dots, 
no dominance of a particular channel occurs. 

Figure 3. Same as in Figure 2, with all sfermion masses fixed at mj = 
3 TeV. 

Figure 4. Scatter plots for neutralino relic abundance Q x h 2 as a 
function of the neutralino mass m x . M 2 and fx are varied in the ranges 
20 GeV < M 2 < 6 TeV and 20 GeV < \fi\ < 3 TeV; tan/3 is varied in the 
range 6 < tan/3 < 10; lightest scalar Higgs boson mass is rrih = 50 GeV; 
sfermion masses are taken as in Figure 2. (a) and (b) refer to neutralino 
compositions which are dominantly higgsino (a\ + a 2 < 0.1) or dominantly 
gaugino (a± + a 2 > 0.9), respectively; (c) refers to the case when higgsino 
and gaugino components are maximally mixed (0.45 < a\ + a 2 < 0.55), for 
positive fx; (d) the same as in (c), for negative fi. 

Figure 5. Same as in Figure 4(a and b), except for sfermion masses 
fixed at mj = 3 TeV. 

Figure 6. Scatter plots for neutralino relic abundance O x /i 2 as a 
function of the neutralino mass m x . M 2 and fx are varied in the ranges 
20 GeV < M 2 < 6 TeV and 20 GeV < \fi\ < 3 TeV; tan/3 is fixed at the 
value tan/3 = 2; lightest scalar Higgs boson mass is rrih = 80 GeV; sfermion 
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masses are taken as in Figure 2. (a) and (b) refer to neutralino composi- 
tions which are dominantly higgsino (a\ + < 0.1) or dominantly gaugino 
(of + ol > 0.9), respectively; (c) refers to the case when higgsino and gaugino 
components are maximally mixed (0.45 < a\ + a\ < 0.55), for positive fx; (d) 
the same as in (c), for negative fi. 

Figure 7. Same as in Figure 6(a and b), except for sfermion masses 
fixed at mj = 3 TeV. 

Figure 8. Same as in Figure 6, except for tan/? = 20 and = 50 GeV. 
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